Competition Among Companies: Coexistence and Extinction 
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We study a spatially homogeneous model of a market where several agents or companies compete 
for a wealth resource. In analogy with ecological systems the simplest case of such models shows 
a kind of "competitive exclusion" principle. However, the inclusion of terms corresponding for 
instance to "company efficiency" or to (ecological) "intracompetition" shows that, if the associated 
parameter overcome certain threshold values, the meaning of "strong" and "weak" companies should 
be redefined. Also, by adequately adjusting such a parameter, a company can induce the "extinction" 
of one or more of its competitors. 

PACS numbers: 89.65.Gh, 87.23.Cc, 05.45.-a 



INTRODUCTION 

During the last few years we have witnessed a wealth of 
work on the application of methods of statistical physics 
to the study of economic problems configuring what some 
authors called econophysics pi 0, @, Within this 
framework a great deal of effort was dedicated to the 
analysis of economic data ranging from stock exchange 
fluctuations 0, production models, size distribution of 
companies, the appearance of money, effects of control 
on the market to market critical properties 0. Another 
problem that attracted enormous interest was the origin 
of power (Paretto) laws, and lognormal distribution with 
power law tails, for the income of individuals, wealth dis- 
tribution, debt of bankrupt companies |(| . An interesting 
source of several mathematical descriptions and models 
used in economical and sociological contexts can be found 
in§. 

Here, our interest is the study, in a deterministic way, 
of aspects of the competition and coexistence of agents 
or companies in a common market. We present a sim- 
ple "toy" model describing, in analogy to some ecological 
problems, a situation of competition among several com- 
panies. Within a "Malthusian-like" model we analyze the 
effect of a kind of intracompetitive contribution on the 
possibility of companies coexistence in a certain market, 
and the changing leadership role (measure through some 
wealth parameter) between, according to the standard 
ecological definition, "strong" and "weak" companies. 

According to ecological studies, starting with 
Volterra's first results on the mathematical theory 
of competition ||, the problem of competition and 
coexistence between species has been analyzed and 
resumed within the Competitive Exclusion Principle (or 
Ecological Theorem), that states: N species that compete 
for n{< N) food resources, cannot coexist |9j. Several 
aspects of this problem have been analyzed by different 
authors, emphasizing, for instance, the conflict between 
the need to forage and the need to avoid competition; 
effects of diffusion-mediated persistence. Generally, 



the system describing competition between species can 
be represented by a set of N differential equations for 
the species and n for the resources. As an example, 
with only one food resource we find only two stable 
stationary solutions: the trivial one (extinction of all he 
species), and that corresponding to the survival of only 
one species, the "strongest" one. There are also studies 
of the problem related to the possibility of coexistence 
in the form of wave-like solutions (ll], . 

Here we adapt the model used in Refs.Jll], [l2| for the 
case of a homogeneous competitive market with a unique 
wealth resource and several firms. In the next Section 
we introduce the model and some particular, instructive, 
solutions. Due to the difficulties of finding analytical 
solutions for the general case, in Section III we focus on 
a representative system with a small number (here 4) of 
firms, and analyze its behavior by numerical methods. In 
the last Section we draw some conclusions. 



THE MODEL 

We start describing the model we use, which is related 
to the one used in Refs.Jll], [l2| for the study of coexis- 
tence in an ecological framework in the form of wave-like 
solutions. Such a model has been adapted to the prob- 
lem of competition of N agents or companies for a unique 
common wealth resource. We indicate with Uj the "size" 
or wealth-parameter representing the welfare of the j- 
company (j — 1, N) and with M the total amount of 
"wealth" . The set of differential equations that we use to 
describe the behavior (for the homogeneous case) of such 
a complex system includes a "Malthusian-like" birth- 
death equation |L3| for each company (f3j M (t) rij (t) cor- 
responding to benefits coming from the wealth's share 



and —ajUj(t) for the "standard" losses-or costs- of the 
j-th company). We also include a contribution that cor- 
responds, in ecological language, to taking into account 
the existence of some kind of intracompetition, that is if 
the j-th the company is alone and can get all the wealth 
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M, it can only grow up to a maximum bounded size, a 
behavior that can be modelled by a Verhulst-like term 
This takes into account the in ecological contexts 
so called carrying capacity of the economic environment 
fhi] , p"5| . In economic terms it could correspond to an in- 
crease of the company efficiency, for instance, through the 
reduction of operative internal costs, improved manage- 
ment, avoiding of competition between different branches 
of the same firm, instrumentation of new technologies, 
etc. 

In addition, instead of assuming a constant finite 
wealth resource as in the so called "Total Wealth Con- 
served model" we consider that M has its own dy- 
namics. For the equation for M, the economic wealth 
accessible to the companies that, in order to simplify 
this initial analysis we assume is unique, we also consider 
a "Malthusian-like" behavior including: its production 
(new resources and technologies, harvest and grain pro- 
duction, etc) that we assume has a constant rate Q, and 
its disappearance due to: (a) natural degradation or rot- 
ting of crops, technologies becoming old, some resources 
being exhausted, that we assume has a rate proportional 
to the total wealth amount —GM (only a certain portion 
of M disappears), (b) the share of each company given 
by -f3 k n k (t)M(t). 

The set of equations is (j — 1, N) 



dt 
d 
dt 



n 3 (t) = [p j M(t)~a J ]n J (t)- lj ^ 

M(t) = Q(t)-[G + J2Pkn k {t)]M(t). (1) 



Similarly to what was discussed in |Q for the case 
of only two species, we can here define a hierarchy from 
the "strongest" (largest ratio between wealth share and 
standard losses, i.e. largest I 1 ) to the "weakest" (small- 
est ratio) companies. Assuming the following hierarchical 
order 
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we have that n\ is the strongest company while njv is the 
weakest one. It is worth noting in passing that Eqs.(Q) 
resembles the form of multimode laser systems p7[ , mak- 
ing it possible to transfer some results from one system 
to the other. 

The stationary solutions results from taking 4fM = 
and f t nj 



M, 



(j = 1,...,N). We found 
Q(t) 
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■1i 



M s 



The last equation implies one of two possibilities 
. „. _ [ft M * 



0, 



ctj] M s 
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(3) 
(4) 

(5) 



It is clear that for large N, to find the solution of this 
system is not easy. In order to fix ideas we consider 
the simplified case where, instead of the above indicated 
hierarchy, we have that all companies are equivalent, that 
is 



0j = 0, a.j 



7j = 7, 



implying 



In this case we have 



Vj. 



(a/Q)(G + NPn s )]Q 2 



7 (G + /V/3n s ) 2 
that can be rewritten as 



" Q 2 Us 
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a(3N -/G 2 

~Q~ + ~Q 2 



(0) 



-05— ) = 0- 



(7) 

It is possible to find under which conditions at least one 
solution of Eq. (0) is real. However, it is more instructive 
to look for the behavior at small n s (n s ~ 0) as all the 
coefficients of n v s with v > are negative. Hence, we can 
easily obtain that a solution, given by 

I3Q 2 - aGQ 



s afiNQ + 7 G 2 ' 

exists (is positive) if (3Q > aG. In this case, as one of the 
associated eigenvalues is zero, a linear stability analysis 
does not give a clear information about the stability of 
the solution and it is necessary to resort to a more refined 
analysis. 

Another instructive case is to consider 



^1 > & 



CtN 



Here we reduce to essentially the same situation studied 
in [p[ [l2|. In particular, it is coincident with the situa- 
tion studied in [T^ |, but now having an "effective" weak 
species given by (N — l)rij, j = 2, N. As in and as 
discussed in detail latter for the case of several firms, it is 
possible to find a coexistence region when 71 overcomes 
some threshold value. In this case a linear stability anal- 
ysis shows a change in the stability of these solutions. 

As a general analytical study of our system even for 
N not too large is far beyond our interest, in the next 
section we focus on a numerical approach for a case with a 
small, however representative, value of N analyzing some 
relevant situations. 



NUMERICAL RESULTS 

As indicated before, here we focus in a case with TV 
small (in fact N = 4) that shows all the relevant aspects 
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FIG. 1: Asymptotic values of m for different values of 7^~ 1 . 
In both cases pi = 5, p2 = 2.5, p3 = 1.5, p4 = 1. Full line: 
Species 1; Dash line: Species 2; Dotted line: Species 3; Dash- 
dotted line: Species 4. In a) 72 = 1, 73 = 0.5 and 74 = 0.1, 
in b) 72 = 1, 73 = 0.1 and 74 = 0.5. 



FIG. 2: As in Figure 1 with pi = 4, p2 — pa = pi = 2 and a) 
71 = .25, 73 = 0.5 and 74 = 0.1, in b) 71 = 1, 73 = 0.5 and 
74 = 0.1. 



we can expect in the large ./V situation. Denning the ratio 
Pi = a~ we consider several situations: 

a) when the different companies are in hierarchical order, 
that is pi > p 2 > P3 > Pi, 

b) when we have pi > p 2 = pz = pi, 

c) when p x = p 2 > pz = p\. 

Among all the possible scenarios we have chosen those 
showing regions of coexistence of all the species, thus in 
case (a) we investigate the situation with fixed jj, with 
j = 2, 3, 4 and varying 71. That allows us to have a con- 
trol parameter, but we recall that the choice is arbitrary. 
In this case the usual strong and weak concepts make us 
consider the species 1 as the strongest and the 4 as the 
weakest. With the inclusion of the new term we find not 
only the possibility of coexistence but also that the orig- 
inal hierarchical order can be permuted several times as 
the parameter values are varied. Thus, we find that the 
ranking of companies suffers many changes, with com- 
panies interchanging roles several times. Some examples 
of this case are shown in Fig.l. Besides these new fea- 
tures we observe the classical extinction of companies as 
predicted by the exclusion theorem. In all the cases the 
extinction occurs by one species at a time. Figure 1 shows 
two typical results for the stationary values reached by rij 
as function of 7j~ . It is apparent that according to the 
different values, the relative status of each company can 
change and even reversed with several crossing among 
them. In this figure 71 varies continuously from 10 to 0. 
The ratios Pi = % are p\ = 5, p 2 — 2.5, P3 = 1.5, P4 = 1. 
In Fig. la we have 72 = 1, 73 = 0.5 and 74 = 0.1, while 
in Fig l.b 72 = 1, 73 = 0.1 and 74 = 0.5. 

For the case (b) above, we considered as our variable 
72. In this case we have one strong species with p\ = 4, 
and three similar species, P2 = Pz = Pi = 2, that could 



coexist if there was not a strongest one. We want to 
note that though p 2 = pz = Pi the same is not true 
for a and (3 values. Once again, we observe coexistence 
between species and a reorganization of the company or 
agents ranking, not accordingly to the original concept 
of strength but depending of the values of ji . The coex- 
istence is achieved within certain parameter region. The 
extinction is gradual and governed mainly by 7 values. 
The results are shown in Fig. 2, where again we depict 
the stationary values of nj as functions of 7" 1 . In Fig 
2. a we have 71 = 0.25, 73 = 0.5, 74 = 0.1, while in Fig 
2.b 71 = 1., 73 = 0.5, 74 = 0.1. It is apparent that the 
most relevant parameter when considering competition 
is the value of 7*. When species are equally strong or 
weak, we observe that different stationary density levels 
are reached according to ji . If ji is of the same order the 
coexistence is granted. On the contrary, a species with 
a high 7, will not survive even if competing with similar 
species. As an additional feature we observe that a usual 
strong species j will not survive even if competing with 
weaker species if 7^ is much higher than that of the other 
species. 

In case (c) we considered two situations: a first one 
varying 71; and a second varying 73. The results are 
shown in Fig. 3, where again we depict the stationary 
values of nj as functions of 7^" 1 and j^ 1 respectively. 
In both cases 3 p\ = p 2 = 4, P3 = p4 = 2, while in 
a) 72 = 0.25, 73 = 0.001, 74 = 1.5 and in b) 7 i = 0.5, 
72 = 1.5, 74 = 0.01 . We confirm that the coexistence 
can be achieved for proper 7 values. At the same time we 
observe that the usual concept of strongest species cannot 
be applied as in both previous cases, species 2, one of the 
strongest if 73 were zero remains as the weakest species 
due to a high 72 value. 
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FIG. 3: As in Figure 1 with pi = p2 — 4, pz — p4 = 2 and 
a)72 = 0.25, 73 = 10" 3 and 74 = 1.5, in b) 71 = 0.5, 72 = 0.5 
and 74 = 0.01 



CONCLUSIONS 

The results shown above indicates how the study of 
simplified models could help in the understanding of 
the role played by 7, an internal company parameter, 
associated to the company efficiency, in situations where 
the complexity of the economic reality makes very hard 
to obtain a complete model. In our case, the model so 
far studied could give some hints on the behavior of 
systems of companies in competition and the possibility 
of coexistence and the way that one company can use 
to "eliminate" the competing ones adopting a policy 
tending to adequately change its 7. Here we have 
analyzed the effect of explicitly including the carrying 
capacity of the environment within our toy model for 
describing the coexistence of species in competition. 
The results put in evidence the role played by the term 
associated to 7 in the possibility of coexistence. We 
recall that in the absence of such term the exclusion 
principle is valid and only the species (one or more) 
with the highest p survive. It is also this term that, 
within certain parameter region, governs the company or 
agents ranking. A model written in the same terms can 
clearly also be applicable to ecological situations. But in 
this case, rather than a deterministic control of 7, some 
fluctuations or cyclic changes in this parameter should be 
considered. This is the subject of a work of us in progress. 
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